Entanglement in an expanding spacetime 
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We show that a dynamical spacetime generates entanglement between modes of a quantum field. 
Conversely, the entanglement encodes information concerning the underlying spacetime structure, 
which hints at the prospect of applications of this observation to cosmology. Here we illustrate this 
point by way of an analytically exactly soluble example, that of a scalar quantum field on a two- 
dimensional asymptotically flat Robertson- Walker expanding spacetime. We explicitly calculate the 
entanglement in the far future, for a quantum field residing in the vacuum state in the distant past. 
In this toy universe, it is possible to fully reconstruct the parameters of the cosmic history from the 
entanglement entropy. 
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The phenomenon of entanglement has attracted much 
attention in recent years. Its central importance in 
the exciting discipline of quantum information science is 
undisputable: it has emerged as a fundamental resource 
in quantum communication quantum cryptography 
2| , quantum teleportation [3( and quantum computation 
^. Recent effort has begun to translate some of the 
aforementioned concepts to the special relativistic set- 
ting [a, Ifil LI 13 ISli LLil and recently progress has been 
made in examining teleportation 111) and entanglement 
between modes of a quantum field [12| when one of the 
observers is uniformly accelerated. 

However, quantum information remains a discipline 
that commonly avoids the conceptual challenges implied 
by one of the most fundamental insights of modern 
physics, namely that spacetime is dynamical and curved. 
When first attempting to understand the basic princi- 
ples of quantum information, the simplifying assumption 
of a flat, or even non-relativistic, spacetime is justifiable. 
However, since our ultimate goal is to properly under- 
stand the nature of the universe on a more complete 
level, and not merely a restricted, simplified aspect of 
it, we must necessarily extend our analysis to the more 
general curved spacetime domain in which we live. More- 
over, while for many physical systems of interest to quan- 
tum information theory a non-relativistic approximation 
is sufficient, the latter entirely fails for massless parti- 
cles, such as photons, or in the presence of strong grav- 
itational fields. The advent of precision measurements 
in cosmology |l3j | , for instance, produces a host of data, 
whose interpretation will be further enriched by trans- 
lation of concepts from quantum information theory to 
curved spaces. 

In this brief note we consider the effect that an ex- 
panding universe has on the entanglement shared be- 
tween scalar particles residing in that spacetime, demon- 
strating that such a dynamic background structure ac- 
tually creates entanglement. Furthermore, we explicitly 
demonstrate the fascinating possibility to deduce cosmo- 



logical parameters of the underlying spacetime from the 
entanglement shared between certain modes of a quan- 
tum field. 

The observation that entanglement is affected by 
the underlying spacetime structure may appear rather 
startling to those unfamiliar with quantum field theory 
on curved spacetime, but is an immediate consequence 
of the latter 0|. Indeed, many of the constructions of 
quantum theory in flat spacetime, such as the notion of a 
particle, only possess limited validity in the general set- 
ting. Despite such conceptual challenges, we show that 
the utility of entanglement can be fruitfully extended be- 
yond its usual domain of non-relativistic quantum infor- 
mation. 

We present the required theory in fair generality, but 
for simplicity discuss as an explicit example a toy uni- 
verse in two spacetime dimensions. More realistic exam- 
ples differ not in principle, but only in their analytical 
complexity. 

On a curved d-dimensional spacetime M with metric 
g, we consider a complex valued scalar field cj>, whose 
dynamics is given by variation of the action 



S = 



d d x^[g^(x)d^d^~m 2 (b 2 ] 1 (1) 
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where g = (— l) d+1 det g^ and m is a real positive pa- 
rameter. The corresponding equation of motion is 



(□ + m 2 )(/>(x) =0, 



(2) 



where D<p — d^^Jgg^" 'd v (j)) / \/g. Due to the linearity in 
(j>, the space of solutions constitutes a vector space, which 
is made into a Hilbert space 7i by equipping it with the 
time-independent inner product 



(3) 



where E is a spacelike Cauchy surface. As in any complex 
Hilbert space, we can find an orthonormal basis of the 
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form {u p ,u*} such that 

(u p , u q ) = -(u*, u*) = S pq and (u p , u*) = 0, (4) 

where p, q are in some (possibly continuous) index set 
V . In general, no particular splitting of the basis into u p 
and u* is distinguished. In a spacetime with time trans- 
lation symmetry, however, we have a timelike Killing vec- 
tor field K, satisfying Ck9 — 0, where C denotes the Lie 
derivative. The corresponding energy conservation allows 
one to meaningfully classify solutions of @ into positive 
and negative frequency solutions, if they are eigenfunc- 
tions of iCk with positive or negative eigenvalues, re- 
spectively. 

Canonical quantization of the theory consists of pro- 
moting the field <f> to an operator field 

4>(x) = / [a~u(p) + a±u*(p)}, (5) 
Jv 

where a~ and ap~ act on the bosonic Fock space 

F = <C®H®{HoH)®{HoHoH)... (6) 
by linear extension of their action on H° n , 

n 

a-\ qi )....\q n ) = 5>i>--<Pki>--k»> erM, 

i=l (') 

a+\ qi )....\ qn ) = \p)\ qi )....\ qn ) eH o(n+1 \ 

where the order of the kets is irrelevant due to the sym- 
metric tensor product o. The Fock space T inherits the 
inner product J3J from 7i, from which it immediately fol- 
lows that a+ = (a~ y , where f denotes the adjoint with 
respect to the inner product on T . We will from now on 
denote a~ simply as a p . From the above construction it 
follows that 

[<4» a i] = S pq and i a P' a <i\ = K' a \\ = °- ( 8 ) 

These commutation relations are often stated as the 
essence of the canonical quantization procedure for 
bosons. If the u v are positive frequency solutions with 
respect to some timelike Killing vector field K, the op- 
erators a p and a p may be meaningfully interpreted as 
creation and annihilation operators for a particle excita- 
tion (of mass m, momentum p, and energy \Jp 2 + rn 2 ), 
if the vacuum state |0) = lc © On © • • • £ T is interpreted 
as the no-particle state. Note that from this definition it 
immediately follows that a p \0) — for all p <G V. 

On a generic curved spacetime, there exists no global 
timelike Killing vector field K, so that in general no 
meaningful particle interpretation can be attached to a 
state of the quantum field. However, if in some region 
of the spacetime under consideration, there exists a lo- 
cal Killing vector field, we will exploit the fact that a 
particle interpretation exists in that region. In the ab- 
sence of a global Killing vector field, the same quantum 



state then generically has particle interpretations vary- 
ing between those regions that feature their own indi- 
vidual local Killing vector fields, and thus different pos- 
itive and negative frequency modes. The corresponding 
change from one set of positive and negative frequency 
modes {u p , u*} to another set {u p , u*} is simply a change 
of the Hilbert space basis. As creation and annihilation 
operators are defined with respect to a specific mode de- 
composition, changing the latter induces the Bogolubov 
transformations 

a p= [a* pq a q - /3*g<4], (9) 

JqEV 

parameterized by the projection coefficients a pq = 
(iip, u q ) and p pq = — (m p , u*) of the basis change in Ti. In 
the case of non- vanishing (3 pq , it follows that the vacua 
|0) and |0), defined with respect to the different mode 
decompositions, are inequivalent. As a consequence, the 
particle concept, so familiar and widely used in discus- 
sions of conventional quantum information, is a more in- 
tricate one on curved spacetime. 

The majority of explicit calculations in quantum field 
theory on curved spacetime are notoriously difficult. 
However, there exist certain specific models of the space- 
time structure that are exactly analytically soluble. We 
shall consider one such simple model that is asymptot- 
ically flat in the distant past and far future, as this al- 
lows us to easily illustrate the main principles and conse- 
quences of examining entanglement in curved spacetime 
without having to resort to approximations or numerical 
solutions. 

Specifically, we consider a two-dimensional Robertson- 
Walker expanding spacetime with line element 

ds 2 = C{t){cLt 2 -dx 2 ), (10) 

where r is the conformal time and the conformal scale 
factor is given as 

C(r) = l + e(l+tanho-r), (11) 

with positive real parameters e and a, controlling the 
total volume and rapidity of the expansion. This de- 
scribes a toy universe undergoing a period of smooth ex- 
pansion. In the distant past and far future, the spacetime 
becomes Minkowskian since C(r) tends to l + 2e and 1 as 
t — > ±oo, respectively. As a consequence, the vector field 
K = d/dr has the Killing property in both the asymp- 
totic in-region (r — > — oo) and out-region (r — * oo), but 
not for finite r. In the asymptotic regions it is possible to 
sensibly discuss the particle content of a scalar field; in 
the intermediate region, however, the concept of a parti- 
cle breaks down. 

In order to find the solutions to the Klein-Gordon equa- 
tion Eq. J2J on this spacetime, we note that C(r) is in- 
dependent of x. We exploit the resulting spatial transla- 
tional invariance and separate the solutions into 

$p{T,x) = {2*y 1 / 2 e^ X p{T), (12) 
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so that V = K in this example. Inserting this into the 
Klein-Gordon equation, we obtain a simple differential 
equation for \ p (r) which may be solved in terms of hyper- 
geometric functions [l^. We then apply the Lie deriva- 
tive Lk to the solutions <p p , in order to identify the nor- 
malized modes u p which behave like positive frequency 
modes in the remote past, and the positive frequency 
modes u p in the far future, respectively: 



u p (t, x) r - 
u p (t, x) t - 



•+oo 



(4 7 rw in )- 1 / 2 e l(pa; - 

(47TW out ) 



1/2 i(px-uiautr) 



(13) 



where the angular frequencies have the form 

^ in = tf + m 2 ] 1 ' 2 , 
^out = b 2 + m 2 (l + 2e)] 1 / 2 , 
1 . 

W± = -(WoutiWin). 



(14) 



We shall henceforth denote all quantities related to the 
in-region using a bar, and those referring to the out- 
region without a bar. A consequence of the linear trans- 
formation properties of hypergeometric functions is that 
the Bogolubov transformations associated with the trans- 
formation from {u p , u*} to {u pi u*} take the simple form 



= a s a s - p s a 



(15) 



so that mixing occurs only between states labelled by s 
and — s. 

Now consider the case where the ingoing scalar field is 
in a vacuum state risGR+ |0} s |0) _ s , with no excitations 
present in any of the modes (from the point of view of an 
inertial observer in the in-region) . Because of the simple 
mixing properties of the Bogolubov transformations, we 
may focus solely on a component |0)fe|0)_fc of the input 
state (disregarding all other modes, by tracing the to- 
tal density matrix over them, yields an overall factor of 
unity), and express this component in terms of the modes 
in the out-region. As a pure state of a bi-partite system, 
this can be written as a Schmidt decomposition 



|0) fc |0)_ fc = 2 , c n \n)k\n)-k, 

71=0 



(16) 



where n labels the number of excitations in the field mode 
k (as seen by an inertial observer in the out-region) and 
the coefficients c„ are real. An explicit expression for 
the Schmidt coefficients c n can be obtained by applying 
Eq. U3 to Eq. 

OO 

= a k \0)k\d)-k = (a* k a k - ftaLk) ^ c n \n) k \n)- k . 

(17) 

A simple relabelling of the mode excitation number then 
allows one to deduce that 



ft , 



(18) 



whilst taking the inner product of Eq. I|16|) with its her- 
mitian conjugate yields the following value for the first 
Schmidt coefficient: 



c = 1/1 



(19) 



Thus Eq. (|l(jf) shows that a state which is interpreted 
as a vacuum in the in-region appears as a state with 
particle excitations in the out-region. We must interpret 
this fact as the creation of particles as a direct result of 
the cosmic expansion. Recall however that in the interim 
region, when our toy universe is undergoing expansion, 
no sensible notion of a particle exists. 

We are now in a position to apply familiar methods 
of quantum information theory to extract information 
about the entanglement of the state Eq. I|16(l . It is a 
simple matter to construct the asymptotic output state 
bipartite density matrix g = |0)_fc|0)fc fc(0|_fc(0|, describ- 
ing excitations in the scalar field modes k and —k. Be- 
cause the Schmidt coefficients in Eq. ()16fl are non-zero, 
the in-vacuum is entangled from the point of view of an 
observer in the out-region. Since the density matrix g 
describes a pure state, the von Neumann entropy S of 
the reduced density matrix 



Qk 



2_j -k{m\g\m) 



(20) 



presents a well-defined measure for this entanglement of 
the modes k with the modes — k. One finds, after some 
algebra, the entanglement 



S = -Ti(g k log 2 g k ) = log. 



7 



7/(7-1) 



1 



where 



7 



Oik 



sinh 2 (7rii;_/(7) 
sinh 2 (7njj+/er) 



(21) 



(22) 



depends on the cosmological parameters e and er, and the 
momentum k of the selected modes. This means that the 
expansion of the universe creates entanglement between 
massive modes of opposite momenta. Modes of a massless 
quantum field do not get entangled, due to the conformal 
flatness of the model studied here. Although this total 
decoupling of massless modes is an artefact of the specific 
example studied here, it illustrates that the massive and 
massless case are generically qualitatively different, and 
non-relativistic intuition fails entirely. Since < 7 < 1, 
the entanglement is monotonically increasing in 7, and 
we may invert Eq. I|21() to obtain ^(S). 

Indeed, the entanglement between the field modes en- 
codes the entire information about the underlying space- 
time. For light particles, a direct relation between the 
cosmological parameters and the degree of entanglement 
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can be obtained. To see this explicitly, assume that a 
universe of the type discussed above, with cosmological 
parameters a and e, is populated by a particle species 
of mass m -C 2cre -1 / 2 . Then we can consider quanta of 
energy E p — y/p 2 + m 2 such that my/e <C E p -C 2er, and 
we obtain the frequencies 



m 

w+TTT^-e <C cr, 



(23) 
(24) 



such that the cosmological expansion parameter e can be 
determined, by expanding Eq. (|22(1 to leading order in e, 
as a monotonically increasing function of the entangle- 
ment S: 



9 P 2 



(25) 



The cosmological parameter a can be determined from 
the respective entanglement of two modes of slightly dif- 
ferent energy. More precisely, one finds by differentiating 
Eq. (|2*2l with respect to the particle energy E that 



1 + 7(5) 



2 \-f^(s)-i 



E. 



(26) 



using the above approximations. Thus, in our simple toy 
universe, the entanglement of massive states carries the 
complete information about the cosmological parameters. 

In a more realistic four-dimensional setting, neutri- 
nos present a natural candidate for such very light 0] 
and approximately only gravitationally interacting par- 
ticles. The required Bogolubov transformations for a 
four-dimensional Friedmann-Robertson- Walker universe 
have been calculated 0, and for spacetime manifolds 
admitting a spin structure, the Dirac equation for spin- 
1/2 fermions can be written down and quantized. Al- 
though this more realistic setting is analytically consid- 
erably more involved, the qualitative features are the 
same. In particular, DeWitt [ljj gives the general vac- 
uum to multi-particle production and annihilation ampli- 
tudes needed in order to generalize our simple example. 
The degree of entanglement of neutrinos therefore likely 
carries information concerning our cosmic history, a fas- 
cinating speculation. 

Our analysis assumed that the quantum field is in a 
vacuum state from the point of view of an inertial ob- 
server in the distant past. While without further assump- 
tions this presents a natural departure point, one could 
instead calculate the fate of a non- vacuum state, or even 
an initially entangled state. The formalism employed 
above is directly applicable also to those cases. We have 
neglected, however, any back reaction of the quantum 
field on the spacetime through the Einstein equations. If 
the dynamics of the universe are not significantly driven 



by the quantum field under consideration, this provides 
a reasonable approximation. 

In conclusion, by extending our understanding of en- 
tanglement to a curved spacetime background, we have 
learned that the latter has a significant effect. We showed 
that an expanding spacetime generates entanglement be- 
tween certain modes of an only gravitationally interacting 
scalar field. While this is not an unexpected result within 
quantum field theory on curved spacetime, we made the 
interesting observation that, conversely, information on 
the underlying spacetime can be recovered from the en- 
tanglement of very light particles. We exemplified this by 
pointing out how all cosmological parameters of a toy ex- 
panding universe can be extracted from quantum correla- 
tions. We will show elsewhere in detail how our treatment 
can be applied to other spacetimes that are either more 
realistic or otherwise of interest. Far from being of mere 
interest to quantum information theory, a further explo- 
ration of entanglement in a curved spacetime therefore 
promises to contribute to a cross-fertilization between 
gravity and quantum information. Recent progress on 
generally covariant quantum field theory supplies valu- 
able new tools for such investigations |lSj . 
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